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摘 要 : 农作物 种 植 结构 是 农作物 空间 格局 的 重要 组 成 部 分 ,是 区 域 土地 资源 和 水 资源 优化 配置 的 基础 。 以 新 疆 


开 一 孔 河 流域 农业 区 为 研究 区 域 ,综合 利用 作物 物候 信息 和 2016 年 的 MODIS NDVI 时 序曲 线 , 获 得 不 同 作物 生长 


差异 明显 的 关键 期 ,选择 关键 期 的 Landsat 8 OLI 影像 ,构建 主要 作物 提取 知识 规则 ,基于 决策 树 方法 开展 农作物 的 
分 类 识别 。 开 一 孔 河 农业 区 2016 年 主要 作物 种 植 面积 为 5.07 x 10° hm ,其 中 棉花 种 植 面积 最 大 ,为 1.97 x 10° 
hm ,玉米 ,小麦 次 之 。 博 斯 腾 湖 和 开 都 河 农业 区 以 辣椒 、 玉 米 和 小 麦 为 主要 作物 ,种 植 结构 比较 零散 ; 孔 余 河 农业 


区 种 植 结构 比较 单一 ,以 棉花 和 香 梨 为 主要 作物 。 与 仅 利 


] 时 间 序 列 的 MODIS 数据 进行 作物 分 类 识别 的 结果 对 


比 表明 ,综合 利用 MODIS 和 Landsat 数据 的 作物 识别 精度 有 显著 提高 ,总 体 分 类 精度 从 62. 58% 提高 到 88.37% , 
kappa 系数 从 0.53 提高 到 0.86。 该 方法 综合 利用 了 MODIS 数据 的 时 序 特征 和 Landsat 数据 较 高 的 空间 分 辩 率 特 


征 , 有 效 地 避免 了 MODIS 数据 空间 分 辨 率 不 足 而 导致 分 类 精度 较 差 的 情况 ,又 避免 了 Landsat 数据 时 间 分 辨 率 不 


足 所 引起 的 时 相 选 择 盲目 性 或 数据 元 余 , 在 干旱 区 农业 种 植 结构 的 提取 领域 具有 一 定 的 应 用 价值 。 
关键 词 : 农作物 ; 种 植 结构 ; MODIS; Landsat 8; 多 源 数 据 集成 ; 开 都 一 孔 人 淮河 流域 


农作物 种 植 结构 描述 的 是 一 个 地 区 或 农业 生产 
单位 内 农作物 的 类 型 及 其 空间 分 布 ,是 农作物 空 
间 格 局 的 重要 组 成 部 分 ,是 区 域 土地 资源 和 水 资源 
优化 配置 的 基础 ” 。 准 确 高 效 地 获取 农作物 种 杆 
结构 和 主要 农作物 产量 信息 仅 对 作物 种 植 结构 的 
优化 和 调整 具有 重要 意义 ” ,对 于 提高 农业 用 水 效 
率 以 及 缓解 干旱 区 水 资源 紧张 也 具有 重要 的 意义 。 

多 种 农作物 通过 一 定 的 种 植 方式 组 合 形成 了 复 
杂 的 农业 土地 系统 ,单一 的 遥感 影像 难以 覆盖 所 有 
作物 光谱 识别 的 最 佳 时 相 , 利 用 时 间 序 列 遥 感 影像 
表征 作物 季 相 及 物候 变化 规律 ,进而 开展 作物 遥感 
提取 的 方法 得 到 了 广泛 应 用 。 高 时 间 分 辩 率 的 
MODIS 数据 在 大 面积 作物 种 植 信息 提取 人 研究 中 的 
MARL ,实现 了 小 麦 . 玉 米 和 大 豆 等 农作物 
的 分 类 识别 和 种 植 结构 的 提取 研究 。 利 用 时 间 序 列 
的 MODIS 影像 进行 作物 信息 提取 能 够 充分 发 挥 高 
时 间 分 辩 率 的 优势 ,但 受 空间 分 辨 率 的 限制 ,只 能 应 


© 


Wis AHA: 2019-05-14; ”修订 日 期 : 2019 -08 -20 


用 于 大 面积 农作物 种 植 区 域 , 在 农作物 种 植 结构 相 
对 复杂 的 区 域 无 法 获得 理想 的 提取 效果 。 中 等 空间 
分 状 率 的 Landsat 数据 在 作物 种 植 信息 提取 的 研究 
eR ae 实现 了 水 稻 等 典型 作物 的 
息 提 取 。 利 用 单一 中 等 空间 分 辨 率 的 Landsat 影 
像 进行 作物 分 类 只 别 的 方法 一 定 程度 上 降低 了 存在 
大 量 混合 像 元 的 现象 ,提高 了 作物 分 类 的 精度 , 却 难 
以 获得 时 间 序 列 较为 完备 和 密集 的 影像 数据 。 
集成 多 时 空 分 辩 率 的 多 源 遥 感 数据 能 够 增强 作 
物 在 物候 和 空间 等 特征 上 的 差异 信息 ,一 定 程度 上 
克服 了 作物 分 类 时 “ 同 物 异 谱 , 异 物 同 谱 ” 现象 的 发 
生 , 有 利于 较 高 精度 的 作物 识别 。 目 前 ,部 分 学 者 利 
用 时 空 融合 技术 构造 高 空间 和 时 间 分 辩 率 的 时 序数 
据 集 ,开展 基于 多 源 遥 感 数据 的 作物 信息 提取 研究 。 
如 陈 梦 露 等 "小 在 我 国 北方 种 植 大 省 河北 、 黑 龙 江 、 
新 疆 典 型 农 区 ,基于 时 空 数 据 融 合 模 型 ,利用 Land- 
sat 和 MODIS 影像 构建 了 高 时 空 分 辩 率 影像 数据 
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集 , 探 讨 时 空 数据 融合 模型 在 农业 信息 提取 领域 的 
适用 性 ;Hao £04] Laura-Piedelobo 等 中 综合 Land- 
sat 和 Sentinel -2 等 数据 实现 了 多 源 数据 的 作物 信 
息 提 取 研 究 。 虽 然 集 成 多 时 空 分辨 率 遥感 数据 在 农 
业 信息 提取 领域 具有 一 定 的 优势 和 适用 性 ,但 目前 
尚 无 一 种 能 够 适用 于 不 同 数据 源 和 应 用 目的 的 时 空 
融合 算法 出 现 ' 。 已 有 的 时 空 融合 算法 都 有 其 各 
自 的 优点 和 局 限 性 ,多 源 遥 感 影像 的 几何 和 辐射 特 
性 差异 的 校正 、 异 质地 表 混 合 像 元 模型 的 非 线 性 散 
射 特征 的 处 理 和 地 物 时 相 变 化 模型 的 建 模 等 过 程 对 
不 同 的 数据 源 和 地 物 的 变化 会 造成 不 同 程 度 的 误 
差 , 而 且 模 型 比较 复杂 ,不 利于 方法 的 推广 应 用 。 

本 文 以 新 疆 开 都 一 孔 和 从 河流 域 (以 下 简称 开 一 
孔 河 流域 ) 农 业 区 为 研究 区 ,综合 利用 高 时 间 分 辩 率 
的 MODIS 影像 数据 获得 作物 差异 明显 的 关键 生长 时 
期 ,选择 同时 期 较 高 空间 分 辨 率 的 Landsat 8 影像 , 开 
展 作物 种 植 结构 的 提取 研究 ,以 期 提供 一 种 简单 有 效 
地 集成 多 源 数据 的 样本 点 选择 和 作物 识别 方法 ,为 干 
早 区 农业 种 植 结构 信息 的 提取 和 应 用 提供 参考 。 


1 研究 区 概况 与 数据 源 
1.1 研究 区 概况 


研究 区 位 于 新 疆 开 一 孔 河 流域 (82°56' ~ 
88°12’E 40°48’ ~43°%20'N) 农 业 区 (图 1) ,主要 包括 


于 斯 腾 湖 农业 区 1 019 km? 、 开 都 河 农 业 区 2 459 
km? 以 及 孔 八 河 农业 区 2 789 km?! 。 开 一 孔 河流 
域 具 有 明显 的 大 陆 性 干旱 气候 特征 ,总 体 特 点 为 日 
照 时 间 长 . 层 夜 温差 大 ,气候 波动 频繁 降雨 稀少 、 获 
发 剧烈 。 农 业 区 降水 主要 集中 在 4 一 7 月 ,最 大 降水 
量 在 20 ~40 mm ,农作物 生长 期 主要 集中 在 3 一 10 
月 。 博 斯 腾 湖 、 开 都 河 以 及 孔 余 河 农 业 区 是 该 地 区 
重要 的 特色 农产品 生产 基地 ,其 中 博 斯 腾 湖 和 开 都 
河 农 业 区 主要 集中 在 和 静 县 和 硕 县 .看 兽 县 和 博 湖 
县 ,主要 种 植 作物 有 小 麦 、 玉 米 .棉花 辣椒 等 ;和 孔 瞧 
河 农业 种 植 区 主要 集中 在 库尔勒 市 和 尉犁 县 ,主要 
种 植 作物 有 棉花 和 香 梨 。 
1.2 数据 源 及 预 处 理 

研究 采用 的 数据 包括 遥感 数据 .主要 农作物 物 
候 信息 和 野外 考察 样本 点 数据 。 

遥感 数据 包括 Landsat 8 OLI 数据 和 MODIS 数 
据 ( 表 1) ,其 中 Landsat 8 OLI 数据 空间 分 辨 率 为 30 
m, 从 地 理 空间 数据 云 平 台 (http://www. gscloud. 
cn/ ) 下 载 2016 年 4 一 9 月 的 7 Æ Landsat 8 OLI 数 
据 ; 从 NASA 官网 (https://ladsweb. nascom. nasa. 
gov/ search) 下 载 的 2016 年 全 年 23 期 分 辨 率 为 250 
m 的 MOD13Q1 以 及 秋 、 夏 季节 的 4 景 500 m 分 辨 
率 的 MOD09Al 产品 数据 。 
Landsat 8 OLI 数据 预 处 理 主要 包括 辐射 定 标 、 
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图 1 研究 区 地 理 位 置 及 样本 点 空间 分 布 示 意图 


Fig.1 Location of study area and spatial distribution of sample points 


R1 遥感 数据 列表 


Tab.1 Remote sensing data list 


传感器 获取 时 间 / 月 -日 
Landsat 8 04 -14 05 - 16 06 -01 07 -27 08 -04 09 -05 09 -21 
MOD13Q1 01 -01 01 -17 02 -02 02 -18 03 -05 03 -21 04 -06 
04 -22 05 -08 05 -24 06 -09 06 -25 07 -11 07 -27 
08 - 12 08 -28 09 -13 09 -29 10 -15 10 -31 11-16 
12 -02 12 -18 
MODO09 Al 06 -09 07-11 08 - 12 09 - 13 
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KARE .影像 裁剪 等 内 容 。MODIHS 数据 已 经 过 严 
格 的 几何 校正 ,辐射 定 标 、 大 气 校正 等 预 处 理 ,利用 
MRT 工具 对 数据 进行 投影 变换 (正弦 投影 转 UTM_ 
Zone_ 45N) .格式 转换 (HDF 转 TIF) 波段 提取 (ND- 
VI 波段 和 对 应 的 QC 波段 ) .影像 裁剪 以 及 S - CUB 
yy 8-1 yh Be 
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从 中 国 气 象 数 据 网 (https://data. cma. cn ) 获 
得 新 疆 开 一 孔 河 流域 小 麦 .玉米 .棉花 3 种 农作物 物 
修 期 资料 ;实地 调查 获得 辣椒 . 香 梨 和 芦苇 物候 期 资 
料 ( 表 2)。 由 于 博 斯 腾 湖 农业 区 分 布 有 大 量 的 芦 
苇 ,也 将 其 作为 农作物 分 类 识别 的 对 象 之 一 。 

研究 区 内 农作物 样本 点 数据 主要 通过 实地 调查 


表 2， 开 一 孔 河流 域 主要 农作物 物候 候 历 
Tab.2 Main crop phenology of the Kai-kong River Basin 


月 2 3 4 5 6 7 8 9 10 

名 | 下 | 十 | 中 下 | 下 | 中 | 下 | 上 | 中 | 下 | 上 | 趾 | 下 I 上 | 中 | 下 | 上 上 | 中 | 下 || 上 | 中 | 下 | 上 下 
小 麦 | 播种 出 苗 aye BOW BOW ATE FLA 收割 
玉米 播种 出 苗 拔节 抽穗 FLA 收获 
棉花 播种 出 苗 KIF Se 花 铃 裂 铃 棉花 采摘 
辣椒 播种 出 苗 定植 开花 结果 收获 
ER 返青 开花 生长 果实 成 熟 
FÆ 发 芽 展 叶 生长 孕穗 抽穗 ”| 开花 成 熟 落叶 


并 结合 Google Earth 进行 采集 。2018 年 7 月 对 研究 
区 内 的 农作物 种 植 情况 采用 CPS 定点 及 拍照 的 方式 
进行 了 实地 调查 ,采集 实地 的 农作物 样本 点 。 此 外 ， 
结合 Coogle Earth 上 影像 的 色彩 以 及 纹理 等 特征 增加 
样本 点 ,并 删除 2016 年 与 2018 年 影像 不 一 致 的 野外 
考察 样本 点 。 采 集 的 样本 点 总 计 368 个 (图 1) ,其 中 
小 麦 67 个 玉米 50 个 .棉花 63 个、 辣椒 76 M ERN 
42 个 沪 革 70 个 ,主要 用 于 精度 验证 。 


2 研究 方法 


为 了 解决 研究 区 实际 样本 点 无 法 获取 的 困难 ， 
本 文 基于 Landsat 8 OLI 和 MODIS NDVI 影像 构造 
作物 样本 点 。 进 而 得 到 不 同 作物 的 NDVI 年 际 变化 
曲线 ,分 析 可 得 作物 生长 明显 差异 的 关键 期 ,选择 关 
键 期 的 Landsat 8 影像 ,通过 构造 作物 识别 知识 规 
则 ,基于 决策 树 分 类 方法 开展 作物 的 识别 与 分 类 ,并 
基于 单一 MODIS 影像 数据 进行 作物 识别 与 分 类 的 
方法 进行 对 比分 析 ,探讨 基于 多 平台 遥感 数据 在 农 


(a) 时 期 


业 种 植 结构 提取 方面 的 适用 性 。 
2.1 选取 作物 样本 点 

实验 中 采用 的 是 2016 年 的 数据 ,由 于 无 法 获得 
研究 区 实际 的 农作物 样本 点 ,本文 基于 Landsat 8 和 
MODIS 影像 ,结合 不 同 作 物 的 物候 信息 构造 和 筛选 
作物 样本 点 。 

首先 ,利用 专家 知识 ,根据 不 同 作物 在 同一 时 期 
的 NDVI 差异 及 同 种 作物 在 不 同时 期 的 NDVI 变化 ， 
选择 不 同月 份 的 Landsat 8 影像 建立 初始 样本 点 。 
具体 方法 为 (图 2) :选择 a 时 期 影像 和 时 期 影像 ， 
其 中 ab 时 期 影像 中 的 每 个 小 矩形 代表 一 部 分 像 
元 ,颜色 代表 像 元 的 NDVI 值 ,NDVI 不 同 则 颜色 不 
同 。 在 影像 ab 中 ,可 以 区 分 出 深 绿色 和 白色 为 不 
同 的 作物 ,对 a b 影像 进行 到 加 , 则 可 以 区 分 出 红 
色 WE . 蓝 色 和 紫色 4 种 不 同 的 作物 。 结 合 农 作物 
物候 信息 ,再 对 4 类 作物 进行 区 分 ,确定 具体 为 哪 一 
类 作物 ,直到 选择 出 所 有 作物 样本 点 。 基 于 该 方法 ， 
构造 的 候选 样本 点 总 计 394 个 , 其 中 小 麦 43 个 、 玉 


jas > 


(b) 时 期 
图 2 候选 样本 点 选择 方法 


Fig.2 Candidate sample point selection method 
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米 61 个 棉花 52 个 .辣椒 68 个 EAL 71 个 芦苇 99 
个 。 完 成 候选 样本 点 构造 后 ,利用 MODIS 数据 对 候 
选 样本 点 进行 筛选 。 计 算 394 个 候选 样本 点 的 全 年 
23 期 MODIS 影像 的 NDVI 值 进行 提取 ,得 到 每 个 候 
选 样本 点 的 NDVI 时 序 变化 曲线 ,剔除 每 一 类 候选 
样本 点 中 NDVI 明显 异常 的 曲线 和 不 符合 农作物 物 
候 规 律 的 曲线 ,得 到 符合 作物 生长 NDVI 曲线 的 样 
本 点 总 计 149 个 ,其 中 小 麦 16 个 、 玉 米 29 个 棉花 
15 个 PRE 26 个 AE 30 个 .芦苇 33 个 。 结 合 高 


分 辩 率 的 Google Earth 影像 对 筛选 得 到 的 样本 点 进 
行 验证 ,并 与 2018 年 野外 考察 的 样本 数据 中 确认 未 
发 生变 化 的 地 物 样 本 点 进行 对 比 , 结 果 显 示 最 终 得 
选 出 的 样本 点 符合 实际 情况 。 
2.2 作物 生长 期 NDVI 曲线 分 析 

对 最 终 确定 的 典型 作物 样本 点 NDVI 曲线 求 平 
均 得 到 小 麦 玉米、 湾 椒 、 棉 花 、 香 梨 和 芦苇 的 NDVI 
年 际 变化 曲线 (图 3)。 

由 农作物 生长 期 NDVI 曲线 可 知 ,小 麦 EAR AM 


09 - 一 一 小 麦 一 一 玉米 一 + 一 棉花 一 上 一 辣椒 ”一 x* 一 香 犁 ”一 x 一 芦苇 


NDVI 
=] 
a 


0.2 
—* 


0.1 * 二 x Set S* SS 


1 17 33 49 65 81 97 113 129 145 161 177 193 209 225 241 257 273 289 305 321 337 353 
MODIS 数 据 获取 日 期 ( DOY ) 


图 3 不 同 作物 及 易 混 淆 地 表 和 覆盖 NDVI 年 内 变化 曲线 


Fig.3 Annual variation curve of NDVI covered by different crops and confusing surface 


芦苇 3 种 作物 的 生长 期 早 于 辣椒 、 棉 花 和 玉米 ,这 与 
作物 物候 历 是 吻合 的 ,说 明 所 选择 的 样本 点 符合 实 
际 情况 。 小 麦 . 香 梨 和 芦苇 3 种 作物 在 3 月 中 旬 左 
右 开 始 生长 ,而 辣椒 、 棉 花 和 玉米 3 种 作物 则 在 4 月 
中 旬 左 右 开始 种 植 。 据 此 ,我 们 可 以 选择 3 月 中 旬 
到 4 月 上 名 的 关键 时 期 影像 将 其 分 为 两 大 类 ,然后 
再 进行 细 分 。 
2.3 基于 决策 树 分 类 的 作物 识别 与 提取 
根据 不 同 作 物 NDVI 年 际 变化 曲线 找到 差异 明 
显 的 关键 时 间 节 点 ,并 选择 该 节点 所 对 应 的 同时 期 
较 高 空间 分 辨 率 的 Landsat 8 影像 通过 构建 知识 规 
则 ,可 以 最 大 程度 地 利用 较 少 影像 实现 NDVI 差异 
明显 的 作物 分 类 ,但 对 于 生长 期 较为 相似 的 作物 , 仅 
基于 NDVI 无 法 实现 分 类 。 因 此 ,结合 作物 在 不 同 
波段 的 光谱 差异 进一步 进行 分 类 。 同 时 ,为 了 避免 
其 他 信息 对 农作物 提取 的 干扰 ,需要 预先 提取 出 水 
体 和 其 他 非 植被 区 域 。 
2.3.1 MNDWI 规 则 利用 MNDWI( Modification of 
normalized difference water index ) [20] 对 水 体 进 行 提 
取 。 计 算 公 式 如 下 : 


Green — SWIR 
MND Green + SWIR oe 


IP: SWIR 表示 短波 红外 波段 的 反射 率 值 ， 
Landsat 8 OLI 第 6 波段 Landsat 7 ETM + 第 5 波段 的 
反射 率 值 ;Green 表示 可 见 光 绿色 波段 的 反射 率 值 ， 
Landsat 8 OLI 第 3 波段 .Landsat 7 ETM + 第 2 波段 的 
反射 率 值 。 为 了 避免 冰雪 融化 的 干扰 ,选择 5 月 16 
日 (第 137 d) 9 H 5 日 (第 249 d) 以 及 9 月 21 日 (第 
265 d)3 个 时 期 的 影像 分 别 代 表 夏 初 秋 初 和 秋末 来 
提取 水 体 信 息 ,具体 规则 为 (MNDWIs。 > 0. 1) and 
(MNDWI,,, >0. 1) and( MNDWI,,,, >0.1)。 

2.3.2 NDVI ALR) 本 文 利 用 NDVI( NDVI-Normal- 
ized Difference Vegetation Index) (21] 对 农作物 进行 提 
取 。 计 算 公 式 如 下 : 


NIR - Red 
NIR + Red 


NDVI = (2) 
式 中 :NIR 表示 近 红 外 波段 的 反射 率 值 , Landsat 8 
OLI 第 5 波段 或 Landsat 7 ETM + 第 4 波段 的 反射 率 
值 ;Red 表示 可 见 光 红色 波段 的 反射 率 值 ,Landsat 8 
OLI 第 4 波段 或 Landsat 7 ETM + 第 3 波段 的 反射 率 
值 。 选 择 5 月 16 日 (第 137 d) 7 H 27 日 (第 209 
d) .9 月 21 日 (第 265 d) 三 个 均 有 植被 覆盖 时 期 的 
影像 来 排除 裸 地 ,具体 规则 为 :(NDVIss >0. 1) or 


ar 
ll 


(NDVI,,; >0. 1) or(NDVI,,,, >0.1). 
2.3.3 农作物 知识 规则 的 建立 小 麦 ERMEE 
3 种 作物 的 生长 期 早 于 辣椒 棉花 和 玉米 ,5 月 中 各， 
小 麦 处 于 孕穗 期 , 香 梨 处 于 叶子 生长 期 ,植被 覆盖 度 
很 高 ;而 玉米 ,棉花 和 辣椒 处 于 出 苗 期 ,植被 覆盖 度 
很 低 ,NDVI 值 远 低 于 0.4。 选 择 5 月 16 日 (第 137 
d) 的 关键 期 影像 将 其 分 为 两 大 类 ,分 类 规则 为 : 
(NDVI ns >0.4)。 接 下 来 在 每 一 类 中 进行 细 分 。 

(1) 小 麦 A PER 

小 麦 与 香 梨 和 芦苇 分 离 的 关键 节点 在 7 月 下 
旬 ,此 时 小 麦 已 经 收割, 在 7 月 27 日 (第 209 d) 的 
NDVI 值 远 低 于 0.3 ,而 此 时 香 梨 和 芦苇 NDVI 值 在 
0.5 以 上 ,通过 设置 规则 (NDVL <0.5) 将 小 麦 与 
香 梨 和 芦苇 进行 分 离 。 由 于 香 犁 和 芦苇 的 NDVI Æ 
长 曲线 相似 ,单一 的 基于 NDVI 值 无 法 将 两 者 很 好 
地 分 离 ,因此 ,结合 光谱 规则 进行 分 离 。 经 过 多 次 实 
验 及 反复 对 比 ,选择 蓝 色 波段 (band blue ) 和 近 红 外 
波段 (band NIR) 的 值 作为 光谱 特征 参数 对 香 梨 和 
芦苇 的 空间 分 布 特征 进行 识别 和 提取 ,具体 规则 为 : 
FÆ (Bpr <2 500 or Bp <600) 。 

(2) 玉米 AEE JHE 

玉米 与 与 其 他 作物 分 离 的 关键 节点 在 9 月 下 
旬 , 此 时 玉米 已 经 收割 ,棉花 和 辣椒 仍然 处 于 生长 
期 。 在 9 月 21 日 (第 265 d) 玉 米 的 NDVI 值 远 低 于 
0.4, 而 此 时 棉花 和 辣椒 的 NDVI 值 在 0.5 以 上 , 通 
过 设置 规则 (NDVIw <0.5 ) 将 玉米 与 棉花 和 辣椒 
进行 分 离 。 由 于 芦苇 生长 状况 参差 不 齐 ,利用 5 月 
16 日 影像 进行 分 类 时 ,玉米 .棉花 和 辣椒 中 也 被 分 
进 部 分 芦苇 。 因 此 ,在 分 离 棉 花 和 辣椒 之 前 需要 把 
芦苇 提取 出 来 。 选 择 5 月 16 日 蓝光 波段 (band 
blue) ,设置 规则 (B,， 友 800 ) 可 将 芦苇 从 棉花 和 车 
PIR 9 月 21 日 ,棉花 处 于 花 铃 收获 期 ,辣椒 
处 于 采摘 期 ,都 已 停止 生长 。 但 由 于 棉花 是 白色 , 辣 
椒 是 红色 ,其 光谱 特征 存在 差异 。 为 了 准确 识别 二 
者 的 空间 分 布 特征 ,选择 蓝光 波段 (band blue ) 和 近 
红外 波段 (band NIR ) 的 值 作为 特征 参数 对 辣椒 和 
棉花 空间 分 布 特征 进行 识别 和 提取 ,具体 规则 为 :六 
PU Bane =550 or Bur <550) o 基于 分 类 规则 构造 的 
决策 树 如 图 4 所 示 。 
2.4 分 类 精度 评价 方法 

混淆 矩阵 是 表示 精度 评价 的 一 种 标准 格式 ,用 
n 行 n 列 的 矩阵 形式 来 表示 。 具 体 的 评价 指标 有 总 
体 精度 OA (overall accuracy) „Kappa 系数 等 。 


总 体 分 类 精度 描述 的 是 对 每 一 个 随机 样本 ,所 
分 类 的 结果 与 参考 数据 所 对 应 区 域 的 实际 类 型 相 一 
致 的 概率 ,计算 公式 如 下 : 


OA = SX/N (3) 


式 中 :N 表示 样本 点 总 数 ;已 表示 被 分 到 正确 类 别 的 
Kappa 系数 是 一 种 定量 评价 遥感 分 类 图 精度 的 
方法 ,计算 公式 如 下 : 


NSR- ENa) 
Kappa = = = 


: (4) 
N -EXX a) 

式 中 :表示 样本 点 总 数 ;X; 表 示 被 分 到 正确 类 别 的 

样本 数 ;X;, 和 针 ,; 分 别 表 示 第 i 行 和 第 i 列 的 总 样 

本 数量 。 


3 结果 评价 与 分 析 
3.1 分 类 结果 


利用 图 4 构建 的 决策 树 , 对 开 一 孔 河 农业 区 的 
主要 农作物 进行 分 类 ,获得 初步 分 类 结果 。 由 于 决 
策 树 分 类 是 基于 像 元 的 作物 分 类 方法 ,其 分 类 结果 
不 可 避免 地 出 现 一 些 像 元 孤岛 等 “椒盐 ”现象 ,利用 
分 类 后 处 理 功能 模块 对 分 类 结果 进行 小 斑 块 去 除 。 
随后 根据 当地 作物 实际 种 植 的 经 验 知识 ,并 借助 
Goolge Earth 影像 ,对 去 除 小 斑 块 后 的 结果 再 进行 适 
当 的 人 机 交互 处 理 。 由 图 5a 可 知 , 博 斯 腾 湖 和 开 都 
河 农业 区 作物 种 植 主要 为 小 麦 . 玉 米 、 湾 椒 及 芦苇 。 
其 中 芦苇 为 自然 生长 ,其 分 布 集中 在 博 斯 腾 湖 西部 
外 围 区 。 而 小 麦 .玉米 .辣椒 分 布 较 零 散 , 没 有 明显 
的 规律 。 相 比 之 下 ,孔雀 河 农业 区 主要 为 棉花 和 香 
A ,种植 比较 集中 ,只 零星 分 布 一 些 玉 米 和 小 麦 。 开 
展 仅 利用 时 间 序 列 的 MODIS 影像 的 农作物 种 植 结 
构 提 取 ,同样 基于 图 5 的 规则 对 阔 值 进行 适当 调整 ， 
经 分 类 后 处 理 得 到 的 农业 种 植 结构 提取 如 图 5b。 
由 于 MODIS 影像 的 空间 分 辨 率 为 250 m ,混合 像 元 较 
多 ,对 比 图 5a 和 图 $b ,直观 上 可 以 看 出 分 类 的 效果 较 
2 ,很 大 一 部 分 农作物 的 像 元 被 分 到 非 植被 区 域 。 
3.2 分 类 精度 评价 

利用 野外 考察 样本 点 ,采用 混淆 矩阵 的 方法 对 
上 述 分 类 结果 进行 精度 验证 ( 表 4) 。 

综合 利用 MODIS 和 Landsat 8 影像 的 分 类 方法 
与 利用 单一 MODIS 影像 的 分 类 方法 相 比 , 除 香 梨 提 
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图 4 农作物 分 类 决策 树 


Fig.4 Decision tree of crop classification 


(a) MODIS 和 Landsat 数 据 分 类 结果 图 F (b) MODIS 数 据 分 类 结果 图 


wn 

作物 分 类 

C | 市 /县 边界 Be 国庆 ee EK BE 辣椒 0 20 40 80 km 
C FERK EE 小 才 (FS Bt 回避 其 他 Ei 


5 基于 决策 树 的 分 类 结果 图 


Fig.5 Classification result graphs based on decision tree 
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表 4 分 类 精度 


Tab.4 Classification accuracy 


/% 


方法 小 麦 玉米 棉花 辣椒 Ges 芦苇 总 体 精度 Kappa 系数 
MODIS + Landsat 78.48 87.50 94.12 85.92 96. 43 97.62 88.37 0. 86 
MODIS 48. 48 38.71 64.71 68. 42 90.91 84.00 62.58 0.53 


取 精 度 略 微 提高 外 ,其 他 作物 的 提取 精度 都 有 显著 
的 提高 。 综 合 利 用 MODIS 和 Landsat 影像 进行 分 类 
的 评价 结果 ,Kappa 系数 为 0. 86 ,总 体 分 类 精度 达到 
88.37% > EPRE AR ,棉花 分 类 精度 较 高 ,分别 
为 97.62% 96.43% 及 94.12% 。 而 小 麦 .辣椒 与 玉 
米 分 类 精度 相对 降低 。 主 要 是 因为 在 30 m 分 辩 率 
的 Landsat 8 影像 上 ,3 种 作物 的 种 植 结 构 零 散 且 复 
RAK ,造成 混合 像 元 的 出 现 。 此 外 ,小麦 分 类 精度 
较 辣 椒 和 玉米 更 低 的 原因 ,是 小 麦收 割 以 后 的 地 块 
上 会 种 植 蕊 莱 等 其 他 作物 ,导致 小 麦收 割 后 的 NDVI 
值 仍 然 处 于 有 植被 状态 ,从 而 产生 误 分 。 而 基于 单 
一 MODIS 影像 进行 分 类 的 评价 结果 Kappa 系数 为 
0.53 ,总 体 分 类 精度 只 有 62.58% 。 其 中 芦苇 和 香 
梨 分 类 精度 相对 较 高 ,分 别 为 84% 和 90.91% ,而 其 
他 作物 的 分 类 精度 则 较 低 , 其 中 玉米 的 分 类 精度 最 
低 , 仅 为 38.71% ,主要 是 大 量 的 玉米 像 元 被 误 分 为 
小 麦 。 由 于 MODIS 影像 的 空间 分 辨 率 较 低 ,在 250 
m 分 辩 率 的 影像 上 存在 着 大 量 的 混合 像 元 ,这 是 导 
致 作物 分 类 精度 低 最 主要 的 原因 。 
3.3 主要 作物 种 植 面积 分 析 

在 综合 利用 MODIS 和 Landsat 影像 进行 分 类 的 
基础 上 ,进一步 对 各 县 市 的 主要 农作物 种 植 面 积 信 
息 进 行 统计 分 析 ( 图 6) 。 


国 小 麦 图 香 梨 四 玉米 口 棉花 图 辣椒 


面积 /103 hm? 


?| 


图 6 各 县 农作物 统计 图 


Fig.6 Crop statistics of each county 


开 一 孔 河 农业 区 主要 作物 种 植 总 面积 为 5. 07 x 
10° hm*。 其 中 棉花 种 植 面积 最 大 ,为 1. 97 x 10° 
hm 。 其 次 是 玉米 和 小 麦 , 种 植 面积 分 别 为 1. 23 x 


10° hm? 和 1.03 x10” hm ,辣椒 种 植 面 积 为 6. 17 x 
10“ hm , 香 梨 种 植 面积 最 小 为 2. 28 x 10* hm 。 由 
图 6 可 知 , 孔 和 淮河 农业 区 以 棉花 种 植 为 主 , 辅 以 香 
梨 、 玉 米 . 小 麦 和 辣椒 等 作物 , 约 占 整个 开 一 了 筷 河 农 
业 区 主要 作物 种 植 面 积 的 43.39% 。 库 尔 勒 市 和 尉 
犁 县 的 棉花 种 植 面积 均 占 绝对 优势 , 香 梨 主要 分 布 
在 库尔勒 市 。 博 斯 腾 湖 和 开 都 河 农 业 区 小 麦 .玉米 
和 辣椒 的 种 植 面积 较为 均衡 , 约 占 整 个 开 一 孔 河 农 
业 区 作物 种 植 面积 的 56. 61% 。 

和 硕 和 静 和 咎 着 地 区 作物 种 植 结构 比较 分 散 ， 
地 块 较为 破碎 ,30 m 分 辨 率 的 Landsat 影像 上 仍然 
存在 一 定 的 混合 像 元 ,农田 地 块 之 间 的 防护 林 像 元 
会 被 误 分 为 农作物 ,使 得 作物 种 植 结构 的 统计 结果 
出 现 一 定 的 偏差 。 小 麦收 割 以 后 的 地 块 上 会 种 植 蔬 
菜 等 其 他 作物 ,导致 小 麦收 割 后 的 NDVI 值 仍然 处 
于 有 植被 状态 ,从 而 产生 误 分 ,使 得 和 硕 Se 
市 小 麦 的 统计 结果 低 于 实际 种 植 信息 。 在 今后 的 工 
作 中 ,可 以 考虑 结合 高 分 辩 率 的 遥感 数据 开展 作物 
种 植 信息 的 提取 。 


4 结论 


(1) 提出 了 综合 利用 作物 物候 期 .多 时 期 Land- 
sat 8 及 时 序 MODIS NDVI 影像 的 信息 样本 点 选取 
方法 ,为 不 易 获 得 样本 点 年 份 和 区 域 农作物 的 提取 
提供 了 较 好 的 解决 方案 。 

(2) 利用 主要 作物 MODIS NDVI 时 序曲 线 , 找 
到 不 同 作 物 差异 明显 的 关键 时 期 ,进而 选择 关键 时 
期 较 高 空间 分 状 率 的 Landsat 8 数据 ,综合 考虑 作物 
的 NDVI 差异 以 及 光谱 差异 规则 ,基于 决策 树 的 分 
类 方法 对 作物 种 植 结构 较为 复杂 区 域 的 作物 识别 能 
够 取得 较为 理想 的 分 类 结果 。 该 方法 可 以 充分 发 挥 
两 种 数据 的 优势 ,弥补 MODIS 数据 空间 分 辨 率 , 及 
Landsat 8 数据 难 易 对 多 种 作物 进行 区 分 的 不 足 , 简 
单 高 效 地 实现 对 作物 种 植 信息 的 提取 。 

(3) 开 一 孔 河 流域 农业 区 主要 作物 种 植 面积 》 
5.07 x10° hm*。 其 中 棉花 种 植 面积 最 大 为 1. 97 x 
10° hm? 、 玉 米 种 植 面积 为 1.23 x 10° hm 、 小 麦 种 植 
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面积 为 1.03 x 10° hm 、 辣 椒 种 植 面 积 为 6. 17 x 10° 
m 、 香 梨 种 植 面积 为 2. 28 x 10* hm 。 博 斯 腾 湖 和 
开 都 河 农业 区 约 占 整个 开 孔 河 农 业 区 作物 种 植 面 积 
的 56.61% ,以 小 麦 、 玉 米 和 辣椒 为 主 , 作 物种 植 较 
为 零散 。 孔 省 河 农业 区 约 占 整个 开 孔 河 农 业 区 作物 
种 植 面 积 的 43.39% ,作物 种 植 比较 单一 ,以 棉花 和 
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Crop planting structure extraction based on remote sensing data in 


Kai-Kong River Basin, Xinjiang 
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Fuzhou 350108 , Fujian , China ; 
2. National & Local Joint Engineering Research Center of Satellite Geospatial Information Technology , Fuzhou 
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Abstract: Crop planting structure plays an important role in crop spatial pattern, and it is the foundation for the op- 
timal allocation of regional land and water resources. In this paper ,taking the Kai-Kong River Basin agricultural area 
in Xinjiang as the study area,a method for crop classification that comprehensively utilizes high temporal resolution 
MODIS images and high spatial resolution Landsat images was proposed. Due to large scope of study area, it is difficult 
to obtain crop sample points uniformly throughout whole research area with traffic and time constraints. MODIS and 
Landsat data for 2016 ,combined with crop phenology data ,were used to construct experimental sample points ;thereby 
providing a better solution for crop extraction in this area wherein sample access is difficult. NDVI time-series curves 
for different crops were constructed based on experimental sample points. Based on the NDVI time-series curves, the 
critical period of crops during growing season were obtained. For these key periods , Landsat 8 OLI images were select- 
ed. Next, extraction knowledge rules for the main crops were constructed ,and identification and classification of crops 
were performed based on decision tree. In 2016, main crop planting area was 5.07 x 10° hm’ of the Kai-Kong River 
Basin agricultural region ,with largest planting area found for cotton (1.97 x 10° hm’) „followed by those for corn and 
wheat. Bosten Lake and Kaidu River agricultural area was dominated by pepper, corn, and wheat, and the planting 
structure was relatively scattered. Planting structure of the Peacock River agricultural area was relatively simple , with 
cotton and pear as main crops. A comparative experiment based on time-series MODIS images for crop recognition and 
classification was also conducted. Results were verified and compared with sample points of field survey. The accuracy 
of crop classification using MODIS and Landsat data was obviously improved as compared with the accuracy of crop 
classification using only time-series MODIS data. Overall classification accuracy increased from 62.58% to 88.37% , 
and kappa coefficient increased from 0.53 to 0.86. The use of high temporal resolution MODIS data and high spatial 
resolution Landsat data can improve the accuracy of crop extraction to a certain extent; this avoided (1) poor classifi- 
cation accuracy caused by the insufficient spatial resolution of MODIS data and (2) phase selection blindness or data 
redundancy caused by the insufficient resolution of Landsat data. Therefore ,this approach has high potential applica- 
tion value in the extraction of crop planting structure in arid areas. 

Key words: crop; planting structure; MODIS; Landsat 8 OLI; Multi-source data integration; Kaidu-Kongque 


River Basin 


